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We present a search for the standard model Higgs boson using events with two oppositely charged 
leptons and large missing transverse energy as expected in iif — >■ WW decays. The events are 
selected from data corresponding to 8.6 fb~^ of integrated luminosity in pp collisions at y/s — 1.96 
TeV collected with the DO detector at the Fermilab Tevatron Collider. No significant excess above 
the standard model background expectation in the Higgs boson mass range this search is sensitive 
to is observed, and upper limits on the Higgs boson production cross section are derived. 

PACS numbers: 14.80.Bn, 13.85.Qk, 13.85.Rm 



INTRODUCTION 

In the standard model (SM), the Higgs boson appears 
during the spontaneous breaking of the electroweak sym- 
metry SU{2) X U{1) that is responsible for the genera- 
tion of the masses of the W and Z bosons. Although 
the SM requires the existence of this neutral scalar par- 
ticle, its mass (Mh) is a free parameter. Direct searches 
at the CERN e+e" collider (LEP) yield a lower limit of 
Mh > 114.4 GeV Q at the 95% C.L. Precision elec- 
troweak data yield, including the latest W boson mass 
requirements from CDF Q and DO Q , constrain the mass 
of a SM Higgs boson to Mh < 152 GeV [1] at 95% C.L. 

In this Article, we present a search for the SM Higgs 
boson in final states containing two oppositely charged 
leptons {££'=efi, ee, or fj.fi, where small contributions 
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from leptonic r decays are also included) and missing 
transverse energy {J^t), using8.6 fb~^ oipp collisions col- 
lected with the DO detector [5| at the Fermilab Tevatron 
Collider. These three leptonic final states are combined 
to produce a result which supersedes our previously pub- 
lished search for Higgs boson production in the oppositely 
charged dilepton and missing transverse energy final state 
based on data corresponding to an integrated luminosity 
of5.4fb~^ 0]. A similar search was published by the 
CDF Collaboration at the Tevatron using 4.8 fb~^ of in- 
tegrated luminosity and by the ATLAS and CMS Col- 
laborations at the CERN Large Hadron Collider (LHC) 
using 4.7 fb~^ and 4.6 fb""'^ of data, respectively [l, 0. 
Using up to 5.4 fb^^ of integrated luminosity, the combi- 
nation of the results from the Tevatron led to the first ex- 
clusion using the H —J- WW decays, excluding the Higgs 
boson beyond the LEP limits, in the mass range from 
162 to 166 GeV at the 95% C.L. Recently, both 

ATLAS and CMS Collaborations have individually com- 
bined all their searches, and the results from ATLAS have 
excluded a Higgs boson in the mass range from 111.4 to 
116.6, 119.4 to 122.1, 129.2 to 541 GeV, while results 
from CMS excluded a Hi ggs boson in the range 127 to 
600 GeV at the 95% C.L. L 
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The primary signal for opposite charge dilepton signa- 



4 



tures with considerable missing energy arises from pro- 
duction of Higgs bosons by gluon fusion gg — H with 
subsequent decay H — >■ WW £^£'1^'. Additional con- 
tributions to this signature come from vector boson fu- 
sion (VBF), qq' — >• qq'H, where the initial state partons 
radiate weak gauge bosons that then fuse to form a Higgs 
boson, and from production in association with a vector 
boson qq' — )- VH = {W/Z)H . The dominant background 
contribution is from diboson production, in particular. 



contributions from non-resonant pp — >■ WW 



processes. Additionally, two types of instrumental back- 
grounds exist: 1) events with mismeasured Ipr in the 
Drell-Yan process pp Zj^* ^ which contribute 

particularly to the ee and /x/j, final states, and 2) events 
with jets misidentified as leptons and photons converting 
to electrons in W boson or multijet production. Although 
such false identification is rare, the resulting backgrounds 
are sizeable as the rates of iy-|-jets and multijet pro- 
duction are significantly higher than that of Higgs boson 
production. Contributions in the /x/i channel from falsely 
identified muons in ly-l-jets events are relatively smaller. 

The following Article first discusses the simulation 
methods used to predict the yields from signal and SM 
background processes. This is then followed by a brief de- 
scription of the DO detector and of the algorithms used to 
reconstruct and identify the objects used in the analysis. 
The event selection and the multivariate techniques used 
to separate the signal from the background are then dis- 
cussed. The different sources of systcmatics uncertainties 
arc then presented, followed by the results of the search 
for the Higgs boson. 



EVENT SIMULATION 



The dominant background processes for the search are 
Z/7*-fjets, VF-|-jets, diboson, tt, and multijet production 
where jets can be misidentified as leptons. Electroweak 
single top quark production is not considered since its 
contribution is negligible. The H^-|-jets and Z-|-jets back- 
grounds are modeled using ALPGEN [2II, with showering 
and hadronization provided by PYTHIA. Diboson pro- 
duction processes {WW, WZ, and ZZ) are simulated 
using PYTHIA. The Z+jets and iy-|-jets processes are 
normalized using the NNLO cross section calculations of 
Ref. [H which uses the NLO CTEQ6.1 PDFs. The Z bo- 
son pt distribution is weighted to match the distribution 
observed in data (25| , taking into account its dependence 
on the number of reconstructed jets. The W boson px 
distribution is corrected to match the measured Z boson 
Pt spectrum [2^ multiplied by the ratio of the W bo- 
son Pt to Z boson pT distributions as predicted in NLO 
QCD [2^. For the search in the ee and e/i channels, 
the ly-l-jets sample includes contributions from events in 
which a jet or a photon is misidentified as an electron. 
For tt production, approximate NNLO cross sections [l^] 
are used, while the NLO production cross section val- 
ues are used for WW, WZ, and ZZ processes [28[. For 



the irreducible background source, WW production, the 
Pt of the diboson system is modeled using the mc@nlo 
simulation [2^ . All MC samples are processed through a 
GEANT simulation of the detector [30|. Recorded detec- 
tor signals from randomly selected beam crossings with 
the same luminosity profile as data are added to the sim- 
ulated detector signals of MC events in order to model 
effects of detector noise and additional pp interactions. 
The simulated background samples are subsequently nor- 
malized to the integrated luminosity. 



Higgs boson signal samples are simulated using the 
PYTHIA [111 Monte Carlo (MC) event generator with the 
CTEQ6L1 parton distribution functions (PDFs) [Hi for 
115 < Mh < 200 GeV in increments of 5 GeV. The 
normalization of these MC samples is obtained using the 
highest-order cross section calculation available for the 
corresponding production process. The cross section for 
the gluon fusion process is calculated at next-to-next- 
to-leading order (NNLO) in quantum chromodynamics 
with soft gluon resummation to next-to-next-to-leading- 
log (NNLL) accuracy For WH, ZH and vector bo- 
son fusion processes, cross section calculations at NNLO 
are used [3, [l3| • AH signal cross sections are computed 
using the MSTW2008 PDF set [I^. The PDF uncer- 
tainties are assessed according to the recommendations 
given in Refs. [l^, [2^. The Higgs boson branching ra- 
tio predictions are from hdecay [2lj. The distribution 
of the transverse momentum (pt) of the Higgs boson in 
the PYTHlA-generated gluon fusion sample is reweighted 
to match the pt as calculated by hqt, at NNLL and 
NNLO accuracy 



DETECTOR AND OBJECT RECONSTRUCTION 

The innermost part of the DO detector Q is composed 
of a central tracking system with a silicon microstrip 
tracker (SMT) and a central fiber tracker (CFT) embed- 
ded within a 2 T solenoidal magnet. The tracking system 
is surrounded by a central preshower detector (CPS) and 
a liquid-argon and uranium calorimeter with electromag- 
netic (EM), fine and coarse hadronic sections. A muon 
spectrometer resides beyond the calorimetry and is made 
of drift tubes, scintillation counters and toroidal mag- 
nets. The DO detector was upgraded in Spr ing 2006 to 
include modifications to the trigger system |3l| as well 
as an additional inner layer of silicon microstrip track- 
in g in stalled near the beam pipe and referred to as Layer 
[32|. The data used for this analysis include 1.1 fb~^ 
collected before these upgrades (Run Ha) and 7.5 fb^^ 
collected afterwards (Run lib). 

Electrons are identified as clusters in the EM calorime- 
ter and are required to spatially match a track recon- 
structed in the central tracking detector. The electron 
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energy is measured from the calorimeter energy deposits 
within a cone of a radius 7^ = y/iArj)^ + {Acj})^ = 0.2, 
in the four layers of the EM calorimeter and the first 
layer of the hadronic calorimeter, where rj and (f> are the 
pseudorapidity [s^ and the azimuthal angle, respectively. 
Due to the different probabilities for jets and photons to 
be misidentified as electrons in the central (CC, \ri\ < 
1.1) and the forward (EC, \ri\ > 1.5) calorimeter regions, 
different criteria are required for the electron clusters. 
The most important of these are: (1) calorimeter isola- 
tion /iso less than 0.15 (CC) or less than 0.1 (EC), with 

/iso = [Etotin < 0.4) - EemCR- < o.2)]/£;EM(7^ < 0.2), 

where E'tot is the total energy in the isolation cone of 
radius TZ = 0.4 and E^m is the EM energy in a cone of 
radius TZ — 0.2; (2) an EM fraction /em larger than 0.9, 
where /em is the fraction of total energy deposited in the 
EM calorimeter, where the energies are measured within 
a cone of radius TZ = 0.2; (3) a track isolation hiso less 
than 3.5 GeV (CC) or less than (-2.5 x \ri\ + 7.0) GeV 
(EC), where hiso is the scalar sum of the pt of all tracks 
originating from the primary pp interaction vertex in an 
annulus 0.05 < TZ < 0.4 around the cluster; (4) a clus- 
ter shape consistent with that of an EM shower; (5) an 
eight-variable electron likelihood Cs that is required to 
be greater than 0.05, where is constructed using the 
variables: /iso, /em, ^iso, the ratio of the electron cluster 
energy to track momentum (E/p), the number of tracks 
within a cone of radius TZ = 0.4, the track-cluster match 
probability computed from the spatial separation and the 
expected resolution, the track distance to the pp inter- 
action vertex at closest approach (dca), and covariance 
matrices that contain variables that relate the energy de- 
positions between various layers of the calorimeter as well 
as the longitudinal and lateral shower development; and 
(6) an artificial neural network trained using information 
from the tracker, calorimeter and CPS detector to further 
reject backgrounds from jets misidentified as electrons. 

Muons are identified by the presence of at least one 
track segment, reconstructed in the muon spectrome- 
ter, that is spatially consistent with a track in the cen- 
tral detector. The momentum and charge are measured 
by the curvature of the central track. The muon can- 
didate must pass quality requirements aimed at reduc- 
ing false matching and background from cosmic rays. 
Muons are required to be isolated. The isolation vari- 
ables are defined to be the scalar sum of the transverse 
energy in the calorimeter (T,E^^°) within an annular cone 
0.1 < TZ < 0.4 and the scalar sum of the transverse mo- 
menta (Epy'*) of tracks within a cone TZ < 0.5 around the 
muon candidate. In the p,^ channel, the isolation vari- 
ables for each muon must satisfy Ep^^ < 0.25 x pl^ and 
j^^caio ^ Q xpl^, pl^ being the momentum of the muon. 
Similarly in the e/i channel, the isolation variables must 
satisfy Ep^'^ < 0.15 x p!^ and E£;f?'° < 0.15 x pi^. For 
and fifi channels, the momentum of the muon track, 
the momentum of the electron track, and the electron en- 



ergy deposit in the calorimeter are not considered when 
calculating isolation variables for the other lepton. This 
prevents the presence of one lepton to spoil the isola- 
tion of the other lepton in events where the separation of 
leptons in (77, </>) space is smaller than TZ = 0.5. 

Jets are reconstructed from energy deposits in the 
calorimeter using Ban iterative midpoint cone algorithm 
with a cone radius of 0.5 Ah jets are required to 

have at least two associated tracks matched to the pp 
interaction vertex. The efficiency of this requirement is 
adjusted in the simulation to match that measured in 
data. Jets can be identified as likely containing b quarks 
(6-tagged) if they pass a selection cut on the output 
of a multivariate (MVA) based b-tagging discriminant, 
trained to separate b jets from light jets 35|. The jet en- 
ergies are calibrated using transverse momentum balance 
in 7-1-jet events [3^. To account for differences in the 
quark/gluon jet composition between the 7-l-jet events 
and the W/Z+^et events, the jet energies are further cor- 
rected in simulated events to match those measured in 
Z-|-jets data. Comparison of ALPGEN with other genera- 
tors [37[ and with the data [38| shows discrepancies in jet 
77 and dijet angular separation. Therefore a data based 
correction allows for a better modeling of these quantities 
in the ALPGEN Z/7*-|-jets samples. 

The IpT is obtained from the vector sum of the trans- 
verse components of energy depositions in the electro- 
magnetic and fine hadronic sections of the calorimeter 
and is corrected for any identified muons. All energy cor- 
rections to leptons and to jets are propagated to the ]^t- 
Data based corrections are applied to MC samples which 
allow for a better modeling of the calorimeter response 
to unclustered objects. 

In order to increase acceptance, all events satisfying 
any trigger requirement from the complete suite of trig- 
gers used for data taking are considered. While most 
of the candidate events in the analysis are selected by 
single- lepton and dilepton triggers, a gain in efficiency 
of up to 20%, depending on the channel is achieved by 
including events which pass lepton-|-jets and lepton+_^T 
triggers. 



EVENT SELECTION 

Candidate events are selected by requiring at least two 
high-pT oppositely charged leptons (e or ^), to originate 
from the same pp interaction vertex along the beam line 
(i.e., within Az=2cm, where Az is the distance between 
lepton tracks along the beam axis, measured at the dis- 
tance of closest approach to this axis). Additional se- 
lections are carried out in two steps, "pre-selection" and 
the "final-selection" . The following section describes the 
selections imposed in each step and also outlines the addi- 
tional kinematic and quality requirements for the search. 
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Preselection 

In the fifi and ee channels, the two highcst-p^ leptons 
are required to satisfy > 15GeV and > lOGeV 
respectively, whereas in the efi channel, p^ > 15 GeV and 
pt^ > 10 GeV are required. Electrons are required to be 
within the acceptances of the EM calorimeter and track- 
ing system (I77I < 1.1 or 1.5 < I77I < 2.5) and muons are 
restricted to the fiducial coverage of the muon system 
I77I < 2.0. In the ee channel, events are rejected when 
both electrons are found in the EC calorimeter as this 
eliminates only a small contribution to the signal which 
has poor signal to background ratio. Additionally, in the 
ee and fifi final states, the dilepton invariant mass Mf-^f,^ 
is required to be greater than 15 GeV. These criteria de- 
fine the "preselection" stage of the analysis. 

To correct for any possible mismodeling of the lepton 
reconstruction and trigger efficiencies, and to reduce the 
impact of the luminosity uncertainty, scale factors are 
applied to the MC samples at the preselection stage to 
match the data. The normalization factors arc deter- 
mined from Drell-Yan dominated samples within a dilep- 
ton mass window of Mi^i.^ € [80, 100] GeV for ee, /i/i and 
Mi^i^ € [57, 75] GeV for e^, and their differences from 
unity are smaller than the luminosity uncertainty. Fig- 
ures [T]-[3] show a comparison between data and the back- 
ground prediction for the distributions of the kinematic 
quantities for each of the dilepton final states after pres- 
election requirements. In the dilepton mass distributions 
shown in Figs. [2^ and the peak in the signal expec- 
tation at Mz originates from ZH associated production 
where the two observed leptons are from the Z boson 
decay. The differences in the widths of the resonance in 
the ee and /iyu channels is due to detector resolution. The 
transverse mass is defined as 

Mt(£,^t) = y'2-/y •^t-[1-cosA(/.(£,^t)], 

and consequently the minimal transverse mass, M™'", 
shown in figs. [2ji and [Sji, is the minimum of the two 
Mt{(-, Ipr) defined for each lepton. 

Jets are considered in this analysis only if they have 
Pt > 20 GeV and \ri\ < 2.4. The preselected samples 
are further subdivided by the number of jets present in 
the event. Dividing the analysis into different jet multi- 
plicity bins significantly increases the sensitivity of this 
search as the signal and background composition change 
between each sample. In particular, 53 — >■ if — >■ WW sig- 
nal processes populate primarily the and 1 jet multiplic- 
ity bins whereas contributions to higher multiplicity bins 
arise mainly from vector boson fusion production and 
associated VH processes which contain additional jets 
in the event. For the background, WW diboson produc- 
tion tends to dominate lower jet multiplicity bins while tt 
events generally contain two jets that are often fo-tagged. 



Subsequent analysis steps are carried out separately for 
events with zero jets, one jet, and two or more jets in 
order to optimally separate signal from backgrounds, re- 
sulting in a total of nine analysis channels {i.e., three 
dilepton final states with three jet multiplicity bins each). 
The jet multiplicity spectrum of the simulated Z/j* sam- 
ple is corrected to match that of the data for each chan- 
nel considered. These corrections are derived within the 
mass windows as described above and have the primary 
effect of improving the ALPGEN modeling of Z/7*-|-jets. 

The number of events for each jet multiplicity bin at 
preselection can be found in Table ID In general, good 
agreement between data and the expected background 
contribution is observed. At this stage, the Z/j* contri- 
bution is the dominant background source. 

Final Selection 

In the ee and /z/x channels, a multivariate discriminant 
is used to remove the dominant Z/j* background present 
in the preselected data sample. The complete details are 
discussed later in this Article. 

As the Z/j* contribution is smaller in the efi channel, 
kinematic selections are instead applied to suppress back- 
grounds after preselection. For the signal, the is not 
aligned with any of the leptons in the final state, while for 
the Z/'y* background processes, the I^t is mostly caused 
by inaccurate measurements of the energies of the lep- 
tons and tends to point in the direction of one of the two 
leptons. Observables that take into account both the ab- 
solute value and the direction of the I^t are M™™ and 
Mt2 , where Mt2 is an extension of the transverse mass 
for final states with two visible and two invisible parti- 
cles [3^. It is obtained as the minimum of the M™'" 
between either lepton and neutrino pair using a mini- 
mization procedure, where the sum of the momenta of 
the neutrinos is varied under the constraint that the sum 
of the momenta of the lepton pair is the missing trans- 
verse energy in the event. The distributions of these two 
observables in the e/i channel after the preselection are 
shown in Fig. 0] for each jet multiplicity bin. The re- 
quirements A'/^"" > 20 GeV and Mt2 > 15 GeV define 
the final selection for this channel. The number of events 
at this selection stage for the e/x state can be found in 
Table m 

INSTRUMENTAL BACKGROUNDS 

The main instrumental background processes for this 
analysis are due to (1) the mismeasurement of I^t in 
Z/7*-|-jets events, (2) the misidentification of associated 
jets or photons in M^-l-7/jets production as leptons, and 
(3) the misidentification of jets in multijet production as 
leptons. 
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FIG. 1: [color online] The (a) dilepton invariant mass, (b) IpT, (c) between the leptons, and (d) minimum transverse mass 
for the e/i channel at the preselection stage. The last bin also includes all events above the upper range of the histogram (a,b,d). 
The signal distribution shown corresponds to a Higgs boson mass of 165 GeV. The hatched bands show the total systematic 
uncertainty on the background prediction. 



TABLE I: Expected and observed numbers of events at preselection in the e^, ee, and /i/_t final states. The signal is for a Higgs 
boson mass of 165 GeV. 
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FIG. 2: [color online] The (a) dilepton mass, (b) Ipr, (c) between the leptons, and (d) minimum transverse mass for 
the ee channel at the preselection stage. The last bin also includes all events above the upper range of the histogram (a,b,d). 
The signal distribution shown corresponds to a Higgs boson mass of 165 GeV. The hatched bands show the total systematic 
uncertainty on the background prediction. 



Z and W Boson Production 

Background contributions from Z bosons are estimated 
using MC simulations. The mismeasurement of Ipx m 
Z/7*+jets events adds a significant source of background 
particularly for tlie ee and /i/i selections, as shown in 
Figs. [33 andEb. 

A W boson decaying leptonically and associated with 
one or more jets or a photon may contribute to the back- 
ground if a jet is misidentified as a lepton or a photon 
overlaps an isolated track or converts into an electron- 
positron pair. The contribution from these backgrounds 
is estimated using MC simulations, and corrections to the 
contributions of jets and photons misidentified as elec- 
trons are derived using data, as explained below. 

An enriched sample of W + 7/jets not overlapping 
with the signal is selected from events passing all the 
selection criteria except that the charges of the two lep- 
tons are required to be identical. This requirement as- 
sumes that the probability of misidentifying a lepton as 



a jet is independent of the lepton charge, and there- 
fore, the like-charge dilepton sample can be used to esti- 
mate background corrections from misidentified leptons 
in the opposite-charge dilepton sample. Corrections are 
obtained separately for initial state radiation jets and 
photons (ISR-y/j) and for final state radiation photons 
(FSR-y) by splitting this control sample into high dilep- 
ton invariant mass {Mg^i^ > 40 GeV) and low dilepton 
invariant mass {Mi-^i^ < 20 GeV) samples where the con- 
tributions of ISR^/j and FSR^ are, respectively, dom- 
inant. These corrections are applied in the ee and e/i 
final states, whereas they are not required in the fj.fi final 
state due to the smaller W + 7/jets contribution. 



Multijet Production 

A high statistics sample of predominantly multijet 
events, where jets are misidentified as leptons, is ob- 
tained from data by inverting certain lepton selection 
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FIG. 3: [color online] The (a) dilepton mass, (b) (c) between the leptons, and (d) minimum transverse mass for the 
fj,fj, channel at the preselection stage. The last bin also includes all events above the upper range of the histogram (a,b,d). 
The signal distribution shown corresponds to a Higgs boson mass of 165 GeV. The hatched bands show the total systematic 
uncertainty on the background prediction. 



criteria. All other preselection criteria are applied in 
order to model the kinematic distributions of the mul- 
tijet background in the signal region. In the ///i channel, 
the opposite-charge requirement for muons is reversed 
and a correction for the presence of non-multijet events 
in the like-charge sample, estimated from simulation, is 
applied. For the e/x and ee channels, the eight-variable 
electron likelihood selection is reversed, and to normal- 
ize the multijet sample to the actual contribution in the 
signal region, the multijet sample is compared to events 
which pass all the signal selections except that a like- 
charge requirement is imposed. This method accounts 
for any kinematic bias introduced from reversing the elec- 
tron likelihood requirement. Since the probability of a 
jet being misidentified as a lepton (Pij) is independent 
of charge, assuming that there is no correlation between 
the charges of the two misidentified leptons in multijet 
events, the like-charge sample has exactly the same nor- 
malization and kinematics as the actual multijet contri- 
bution. Pij depends on the jet multiplicity, and therefore 



the multijet background is estimated separately for each 
jet multiplicity bin. The analysis further assumes con- 
tributions of non-multijet processes are negligible in the 
reversed lepton quality sample. 



MULTIVARIATE ANALYSIS 

A multivariate technique is used to characterize events 
as originating from a Higgs boson signal or from back- 
ground processes and to achieve maximum separation be- 
tween them. A random forest of boosted decision trees 
(BDTs) j4ol | is used to construct a discriminant from kine- 
matic variables, taking into account their correlations. 
The decision trees are trained separately in each of the 
nine analysis channels and for each Higgs boson mass 
hypothesis. To increase the statistics of the available 
simulated signal events, signal samples for neighboring 
mass hypotheses are used for the training of the multi- 
variate discriminant. For example, the training of the 
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FIG. 4: [color online] A/™'" distribution for the e/j. channel in the (a) 0-jet bin, (b) 1-jet bin, and (c) > 2-jet bin. Mt2 
distribution for the e/j. channel in the (d) 0-jet bin, (e) 1-jet bin, and (f)> 2-jet bin. The last bin also includes all events above 
the upper range of the histogram. The signal distribution shown corresponds to a Higgs boson mass of 165 GeV. The hatched 
bands show the total systematic uncertainty on the background prediction 



discriminant for the 165 GeV mass hypothesis uses sig- 
nal samples corresponding to a Higgs boson mass of 160, 
165, and 170 GeV. 



Multivariate Discriminant against Z/f* 

A BDT discriminant is used in the ee and fj,fi final 
states to reject the large Z/'j* background while retain- 
ing a high signal efficiency. This random forest of BDTs 
wiU be referred to as DY-BDT. The DY-BDT is trained 
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for each Higgs boson mass hypothesis and jet multiphcity 
bin, separately for the ee and /i^ final states, to differ- 
entiate between the Zj^* background and all considered 
SM Higgs boson signal events. 

The following input variables are used for the DY- 
BDT: 

(i) lepton 

(ii) invariant mass of the leptons, Mg^i^ 

(iii) azimuthal opening angle between the two leptons, 

(iv) separation in 77, (/> space between the two leptons, 

(v) minimal transverse mass, Af ™™ 

(vi) extended transverse mass, Mti 

(vii) missing transverse energy, 

(viii) smallest and largest of the azimuthal angles, 
between the Ipr and either lepton 

(ix) transverse mass of the Ipj- and the dilepton pair, 

A/t(^1^2,^t) 

(x) special missing transverse energy, Ip^^^'^^^^ , defined 
for object C, which corresponds to either the nearest 
lepton or jet in the event relative to the direction of 
the ^t: 



, special 
T 



X sin[A(^(^T, C)]: otherwise 



(xi) jet pT 

(xii) scaled missing transverse energy defined as 



I scaled 
T 



Y,.^,^ [AEi-^ ■ sin 0i-^ ■ cos A(b (jet, ^t) 



where AE^°^ is a measure of jet energy resolu- 
tion and is proportional to v£?^; the fluctuation 
in the measurement of jet energy in the trans- 
verse plane can be approximated by the quantity 
A£;j<=' • sin6lj«* Q 

(xiii) azimuthal angle between the J^t and the jets, 
A0(^T,jet) 

(xiv) absolute value of the pseudorapidity difference be- 
tween the jets, |A77(ji, j2)|, where ji and j2 are the 
two highest-pT jets in the event 



Variables (i) and (ii) exploit the di-lepton kinematics 
of the event. Variables (iii) and (iv) are related to the 
opening angle between the two leptons and provide dis- 
crimination against SM backgrounds which tend to ex- 
hibit back-to-back topologies. This is not the case for 
Higgs boson decays because of the spin correlation in the 
scalar decay where leptons tend to be aligned in the same 
direction. 

The ^T-related variables (v)-(ix) help distinguish gen- 
uine IpT in the Higgs boson signal from mismeasured I^t 
in Z/^* events. Variable (x) helps to further suppress 
Z/j* events, which populate lower values of ^ fecial 
where a mismeasured lepton or jet tends to align with 
the IpT direction Q. Variables (xi)-(xv) are used in the 
1-jet and 2-jet bins, as appropriate. Since the events are 
categorized in terms of jet multiplicities, variables (xii)- 
(xv) exploit the jet kinematics in the event. 

To reject most of the Z/j* background after the pre- 
selection, events are required to appear in the signal- 
like region of the DY-BDT discriminant. This defines 
the final selection of the ee and fifi final states. The 
threshold varies for each Higgs boson mass hypothesis 
in each jet multiplicity bin and yields a Z/j* rejection 
factor of 0(10-5), 0(10-^), 0(10^2) for the 0-jet, 1- 
jct, 2-jet bins, respectively for all dilepton channels and 
Higgs boson masses. The thresholds are chosen to ob- 
tain similar rejection factors of background events as 
the cut-based analysis employed in the previous publi- 
cation 0. The DY-BDT discriminants for a Higgs bo- 
son mass of 165 GeV are shown in Fig. [5l This figure 
demonstrates that a good separation is achieved between 
the Z/j* background and the majority of signal. How- 
ever it can be noticed that some signal events cannot be 
distinguished from the background and have a very low 
DY-BDT discriminant value. This is primarily due to 
some of the Higgs decay modes which have a signature 
similar to Z/^* background. The numbers of events at 
the final selection stage for the ee and final states are 
shown in Table [Hi 



Final Multivariate Discriminant 

In the final selection step, the signal is separated from 
the remaining backgrounds using an additional random 
forest of BDTs. This final random forest of BDTs re- 
ferred to as FD-BDT, is trained for each Higgs boson 
mass hypothesis and jet multiplicity bin, separately for 
the three dilepton final states using signal and back- 
ground samples, which satisfy the final selection crite- 
ria, to differentiate between all Higgs boson production 
processes and backgrounds. These decision trees use as 
inputs all the variables from the DY-BDT listed above 
with the addition of the following variables: 



(xv) invariant mass of the two jets, Af(ji, J2). 



(i) electron quality likelihood output, £§; for the di- 
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FIG. 5: DY-BDT discriminant for the 0-jet (top row), 1-jet (middle row) and > 2-jet (bottom row) for the ee [left (a,c,e)] 
and fifi [right (b,d,f)] final states. The discriminant shown is trained for a Higgs boson mass of 165 GeV. A final selection 
requirement is applied in the above distributions of 0.35, -0.6, and -0.85 for the ee final state and 0.9, 0., and -0.7 for the fifi 
final state, in the 0-jet, 1-jet, and > 2-jet bins, respectively. The hatched bands show the total systematic uncertainty on the 
background prediction. 



electron channel the lower value of the two electron 
quality likelihood outputs is used 

(ii) a quality criterion based on the number of hits in the 
muon spectrometer characterized in four distinct 
categories; this parameter is referred to as "muon 



quality" and for the dimuon channel the lower qual- 
ity of the two muons is used 

(iii) number of Layer hits in the SMT matched to each 
electron 
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(iv) track isolation variable of each niuon 

(v) the product of charge and pseudorapidity, for both 
leptons £i and £2 

(vi) 6-tag output: the output of a multivariate discrim- 
inant to separate jets originating from heavy flavor 
quarks (6 and c) from those originating from light 



partons; for the channels with iVjots > 2, the small- 
est and largest 6-tag outputs are used. 

Some representative input distributions to the FD- 
BDT at the final selection stage with all jet multiplicity 
bins added in each distribution are shown in Figs. Inland 

m 
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FIG. 7: [color online] The Ip y^'^'^' for the ee (a) and /i/i (b) channel at the final selection stage. The last bin includes all events 
above the upper range of the histogram. The signal distribution shown corresponds to a Higgs boson mass of 165 GeV. The last 
bin also includes all events above the upper range of the histogram. The hatched bands show the total systematic uncertainty 
on the background prediction. 
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FIG. 8: [color online] The quality variable for (a) the ee channel and (b) the /x/i channel at the final selection stage. The signal 
distribution shown corresponds to a Higgs boson mass of 165 GeV. The hatched bands show the total systematic uncertainty 
on the background prediction. 



Representative distributions of the electron and muon 
quality variables, £§ a-nd "muon quality," are shown in 
Fig. [51 These along with other variables given in (iii) 
and (iv) gauge the quality of the reconstruction of the 
lepton and are crucial to discriminate between true lep- 
tons and jets misidcntificd as leptons originating from 
backgrounds like VF+jets. The distribution for the prod- 
uct of charge and pseudorapidity, is symmetric in rj for 
the signal, however this is not true for the background 
processes with misidentified leptons. 

The output from ^-tagging is used to separate the 
Higgs boson signal from ti production, which is an im- 
portant background in the 1 and 2 jet multiplicity bins. 
An MVA-bascd ^-tagging [35| is employed in each of the 
dilcpton final states to discriminate the signal, which 



comprises primarily light flavor quarks, against the heavy 
flavor jets arising from top quark decays. The distri- 
butions for smallest and largest 6-tagging output in the 
> 2-jet multiplicity bin arc shown in Fig. |9l 

The distributions of the final BDT discriminant for 
each channel and Higgs boson masses of 125 GeV and 
165 GeV are shown in Figs. [TOl-IIl 



SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties are characterized for each fi- 
nal state, background, and signal processes. Uncertain- 
ties that modify only the normalization and uncertainties 
that change the shape of the final discriminant distri- 
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FIG. 9: [color online] The (a) smallest 6-tagging MVA output and (b) largest 6-tagging MVA output for the > 2-jet multiplicity 
bin, for all channels summed up at the final selection stage. The signal distribution shown corresponds to a Higgs boson mass 
of 165 GeV. The hatched bands show the total systematic uncertainty on the background prediction. 



TABLE II: Expected and observed numbers of events after the final selection in the e/i, ee, and /i/i final states. The signal is 
for a Higgs boson mass of 165 GeV. 





Data 


Total background 


Signal 


Z/l* 


tt 


W + 7/jcts 


Dibosons 


Multijet 




1729 


1806 


30 


94 


335 


766 


584 


29 


jets 


1117 


1222 


18 


70 


11 


641 


486 


15 


1 jet 


335 


307 


8 


19 


98 


94 


87 


10 


> 2 jets 


277 


277 


4 


5 


226 


31 


11 


4 


ee: 


1607 


1644 


14 


466 


200 


658 


288 


33 


jets 


812 


881 


8 


135 


6 


499 


222 


20 


1 jet 


430 


408 


4 


181 


54 


114 


52 


7 


> 2jets 


365 


355 


2 


150 


140 


45 


14 


6 


fi^: 


1950 


1997 


18 


1101 


231 


198 


328 


140 


jets 


645 


720 


10 


227 


4 


155 


236 


98 


1 jet 


581 


564 


5 


376 


56 


35 


68 


30 


> 2jets 


724 


713 


3 


498 


171 


8 


24 


12 



button are taken into account. Systematic uncertainties 
that contribute only to the normalization are: theoreti- 
cal cross sections for diboson, 6%, and production, 7%; 
multijet normalization, 30%; overall normalization, 4%, 
which accounts for the uncertainty on the lepton trig- 
ger/identication efficiency and the integrated luminosity; 
and a Z+jcts jct-bin-dependent normalization (2-15)%. 

Since the analysis is split into categories depending on 
the number of reconstructed jets, renormalization and 
factorization scale uncertainties on <j{gg — )■ H) are esti- 



mated following the prescription described in Ref. [41 



By propagating the uncorrelated uncertainties of the 
NNLL inclusive [Illii], NLO > 1 jet [l^, and NLO > 2 
jets [ill cross sections to the exclusive gg — >■ H+Q jet, > 1 
jet, and > 2 jets rates, the uncertainty matrix shown in 
Tablelmlis built. The PDF uncertainties for aiog ^ H), 
obtained using the prescription in Refs. [isl . |20| are also 
summarized in Table lllll The uncertainties on the inclu- 



sive a{VH) and a{qqH) are taken as 6% and 5%, respec- 
tively. 

Sources of systematic uncertainty that affect both the 
normalization and the shape of the final discriminant dis- 
tribution are: jet energy scale (1-4)% and jet energy res- 
olution (1-3)%, determined by varying the parameters 
of the energy scale correction and the energy resolution 
function within one standard deviation (s.d.) of the un- 
certainty and repeating the analysis using the kinemat- 
ics of the modified jets; jet association to the pp inter- 
action vertex (1-2)%, obtained by varying the correc- 
tion factor within its uncertainty; shape of the 6-tagging 
discriminant associated with heavy flavor jets (3-5)%, 
determined by varying the correction factor of the b- 
tagging neural network output within its uncertainty; 
M^-l-jets modeling (6-50)%, depending on jet multiplic- 
ity bin and final state, obtained by varying the correction 
factors described in the "Instrumental Backgrounds" sec- 
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FIG. 10: [color online] Final BDT discriminant for the (top-row) 0-jet, (middle-row) 1-jet, and (bottom-row) > 2-jet bins for 
the ejj, final state for a Higgs boson masses of 125 GeV [left (a,c,e)] and 165 GeV [right (b,d,f)]. The hatched bands show the 
total systematic uncertainty on the background prediction. 



tioii within their uncertainties. These uncertainties are 
presented in terms of the average fractional change across 
bins of the final discriminant distribution for all back- 
grounds and depend on the jet multiplicity. 

Several systematic uncertainties are also included 
which have a small (< 1%) effect on the background 
model: modeling of diboson production in terms of 



prp(WW), determined by taking the fractional differ- 
ence of the predicted final discriminant shape between 
MC@NLO and PYTHIA generators; modeling of diboson 
production in terms of the impact of the gluon fusion 
production process on the A(p distribution between the 
leptons, determined by taking the fractional difference of 
the predicted final discriminant shape between MC@NLO 
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FIG. 11: [color online] Final BDT discriminant for the (top-row) O-jet, (middle-row) 1-jet, and (bottom-row) > 2-jet bins for 
the ee final state for a Higgs boson masses of 125 GeV [left (a,c,e)] and 165 GeV [right (b,d,f)]. The hatched bands show the 
total systematic uncertainty on the background prediction. 



and GG2ww [i^l generators; and the pr of the vector 
boson from M^-fjets and Z+jets production. A summary 
of the dominant systematic uncertainties is given in Ta- 
ble |TVl 



RESULTS 

The methodology of this search is validated by an inde- 
pendent measurement of the pp — >■ W^W~ cross section 
using the analysis procedure described in the "Event Se- 
lection" section of this Article, considering WW events 
as the signal. This is motivated by the fact that WW 
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FIG. 12: [color online] Final BDT discriminant for the (top-row) 0-jet, (middle-row) 1-jet, and (bottom-row) > 2-jet bins for 
the final state for a Higgs boson masses of 125 GeV [left (a,c,e)] and 165 GeV [right (b,d,f)]. The hatched bands show the 
total systematic uncertainty on the background prediction. 



TABLE III: Elements of the uncertainty matrix of the scale {^K,^^F) and PDF uncertainties on a-{gg — >■ H) for the three jet 
multiplicity categories considered, where so, si and S2 are the elements of the uncertainty matrix. 





So 




32 


PDF 


Ojet 


13.4% 


-23.0% 




7.6% 


> 1 jet 




35.0% 


-12.7% 


13.8,% 


> 2 jets 






33.0% 


29.7% 
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TABLE IV: Summary of systematic micertainties (in %) for source categories. The jet, 6-tagging and PDF related uncertainties 
are quoted for all the backgrounds combined. 



Source 


Uncertainty (%) 


Overall normalization 


4.0 


W+jets normalization 


6.0-50.0 


Diboson cross section 


6.0 


tt cross section 


7.0 


Multijet normalization 


30.0 


Z+jets jet-bin normalization 


2.0-15.0 


gg ^ H cross section 


See Tablelm] 


VH cross section 


D.U 


qqH cross section 


5.0 


Jet energy scale 


1.0-4.0 


Jet resolution 


1.0-3.0 


Jet primary vertex association 


1.0-2.0 


6-tagging discriminant 


1.0-2.0 


PDF (background) 


2.5 



production is the main contributor to the diboson entry 
in both Tables H] and HI] compared to the expected yields 
from WZ and ZZ production backgrounds. Similarly to 
the Higgs boson search, a dedicated BDT is constructed, 
but now it is trained to separate WW production sig- 
nal from other SM processes. For this BDT, we use the 
identical input variables, the same separation method in 
terms of jet multiplicity bins, and the same treatment 
of systematic uncertainties as in the Higgs boson search. 
The ee and fifi final states use only the and 1 jet mul- 
tiplicity bins while the efi final state uses all three jet 
multiplicity bins yielding a total of seven analysis chan- 
nels for the combination. The results obtained for the 
WW cross section in the individual final states and their 
combination are summarized in Tabic |Vl The measured 
value of 11.1 ± 0.8 pb is in good agreement with the SM 
prediction of 11.7 ± 0.8 pb [28[. The presence of a Higgs 
boson signal in the mass range 115 < Mjj < 180 GeV 
would bias the cross section measurement result by 5% at 
most. This maximum bias is reached for AIh = 165 GeV, 
but at low masses {Mh < 130 GeV), the bias would be 
less than 2%. 

Figures [T51 and [Til show the expected WW and Higgs 
boson signals, respectively, for the combined decay chan- 
nels in the analysis. In these distributions, the data is 
shown, ordered in bins of increasing values of the s/b ra- 
tio, after the subtraction of the SM backgrounds. The 
background model is fit to the data, and the uncertain- 
ties on the background are those after the systematic 
uncertainties have been constrained by the fit. 

The final multivariate discriminants of the SM Higgs 
boson search, shown in Figs. [10] - [l2l demonstrate that 
the data is well described by the sum of the background 
predictions. In the absence of an excess in the num- 
ber of observed events above the SM backgrounds, these 



BDT output distributions are used to set upper limits 
on the Higgs boson inclusive production cross section 
a{pp —T' H + X) assuming SM values for the branching 
ratios and for the relative cross sections of the various 
Higgs production mechanisms considered. The limits are 
calculated using a modified frequentist method with a 
log-likelihood ratio (LLR) test statistic [i^. The value 
of CLg is defined as CLg = CLg+b/CLb, where CLg+t 
and CLf, are the p- values for the signal-|-background and 
background-only hypotheses, respectively. Expected lim- 
its are calculated from the background-only LLR distri- 
bution whereas the observed limits are quoted with re- 
spect to the LLR values measured in data. They both 
are reported at the 95% C.L. 

The multivariate discriminants corresponding to the 
nine individual channels are all used to obtain upper lim- 
its on the Higgs boson production cross section. Given 
the differences in the background contributions to each 
of the channels, the nine BDT output distributions are 
not combined in a single distribution for the limit ex- 
traction, but treated separately. The degrading effects 
of systematic uncertainties on the search sensitivity are 
minimized by fitting individual background contributions 
to the data by maximizing a profile likelihood function for 
the background-only and signal-|-background hypotheses 
separately, taking into account appropriately all correla- 
tions between the systematic uncertainties |46|. Table [VTl 
and Fig. [TS] present expected and observed upper limits 
at the 95% C.L. for a{pp H + X) relative to SM pre- 
dictions for each Higgs boson mass considered. 

The corresponding LLR distributions are shown in Fig. 
[T6l Included in this plot are the median of the LLR dis- 
tributions for the background-only hypothesis (LLRb), 
the signal-plus-background hypothesis [LLRs+b), and 
the observed value for the data (LLRohs)- The shaded 
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FIG. 13: [color online] The background-subtracted data distribution for the final discriminant, summed in bins with similar 
signal to background ratio, for the WW cross section measurement. The uncertainties shown on the background-subtracted data 
points are the square roots of the post-fit background predictions in each bin, representing the expected statistical uncertainty 
on the data. Also shown is the ±1 standard deviation (s.d.) band on the total background after fitting. 



TABLE V: Summary of the measurements of the pp — >■ W^W cross section measurement (in pb) in all seven channels 
considered and their combination. 



Channel 


a{pp^W+W-) 




10.6 ± 0.6 (stat) ± 0.6 (syst) 


ee 


12.4 ± 1.2 (stat) ± 0.9 (syst) 


MM 


11.0 ± 0.9 (stat) ± 0.7 (syst) 


Combined 


11.1 ± 0.5 (stat) ± 0.6 (syst) 



bands represent one and two s.d. departures for LLRh 
centered on the median. The separation between the 
LLRh and LLi?s+6 distributions provides a measure of 
the discriminating power of the search. The current 
result indicates that the signal-|-background model can 
be separated from the background-only model by up 
to 1 s.d. over most Higgs boson masses between 115 to 
200 GeV while the level of separation increases above 2 
s.d. for Higgs boson masses between 160 to 170 GeV. The 
sensitivity of the search reaches an expected exclusion of 
159 < Mh < 169 GeV at 95% C.L. However due to a 
slight excess in the data, an observed exclusion is not 
obtained. 



The best observed limit is obtained at 160 GeV, where 
it reaches 1.1 times the SM expectation. This channel 
is the single most sensitive channel when the H — > WW 
branching ratio is dominant {Mh > 135 GeV), and for 
lower masses at Mh = 125 GeV, this search still has 
a similar sensitivity as a single major low mass channel 
(WH or ZHVwith an expected limit of 3.8 times the SM 
expectation [47[ . The results and the analysis techniques 
arc validated through an independent measurement of 
the WW production cross section, which agrees with the 
NNLO calculation. 



ACKNOWLEDGMENTS 



CONCLUSIONS 

We have performed a search for SM Higgs boson pro- 
duction using final states with two oppositely charged 
leptons and large missing transverse energy in the e//, 
ee, and p,p, channels. After imposing all selection crite- 
ria, no significant excess in data over expected SM back- 
grounds is observed. We set upper limits on Higgs boson 
production at the 95% C.L. The sensitivity of the search 
reaches an expected exclusion of 159 < Mh < 169 GeV. 



We thank the staffs at Fermilab and collaborating 
institutions, and acknowledge support from the DOE 
and NSF (USA); CEA and CNRS/IN2P3 (France); 
MON, Rosatom and RFBR (Russia); CNPq, FAPERJ, 
FAPESP and FUNDUNESP (Brazil); DAE and DST (In- 
dia); Colciencias (Colombia); CONACyT (Mexico); NRF 
(Korea); FOM (The Netherlands); STFC and the Royal 
Society (United Kingdom); MSMT and GACR (Czech 
Republic); BMBF and DFG (Germany); SFI (Ireland); 
The Swedish Research Council (Sweden); and CAS and 



21 



in 


80 


o 
d 


60 


w 40 


c 

0) 


20 


> 







-20 




-40 




-60 




-80 



D0, 8i fb \ W + ^j 




R1 
o 

w 
c 

> 

LU 



■1 



OOM50, 8.6 fb \ II + 
Tb) 



50 

50 
00 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
Final Discriminant 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
Final Discriminant 



FIG. 14: [color online] The background-subtracted data distributions for the final discriminants, summed in bins with similar 
signal to background ratio, for (a) Mh ~ 125 GeV and (b) Mr = 165 GeV. The uncertainties shown on the background- 
subtracted data points are the square roots of the post-fit background predictions in each bin, representing the expected 
statistical uncertainty on the data. Also shown is the ±1 standard deviation (s.d.) band on the total background after fitting. 



TABLE VI: Expected and observed upper limits at the 95% C.L. for o-(pp — s- H + X) relative to the SM for the total combination 
and separately for the e/i, ee and /i/i channels for different Higgs boson masses (Mh)- 



Mh (GcV) 


115 


120 


125 


130 


135 


140 


145 


150 


155 


160 


165 


170 


175 


180 


185 


190 


195 


200 


Exp. 
Obs. 


all: 
all: 


8.00 
13.27 


5.37 
9.14 


3.81 
5.00 


3.02 
4.71 


2.43 
3.93 


2.09 
3.28 


1.77 
2.13 


1.53 
1.99 


1.28 
1.75 


0.92 
1.10 


0.85 
1.17 


1.05 
1.40 


1.27 
1.40 


1.49 
1.64 


1.88 
1.91 


2.48 
2.34 


2.87 
2.87 


3.32 
3.50 


Exp. 
Obs. 


e/i 
efj, 


11.25 
13.86 


7.08 
8.50 


5.07 
5.12 


4.01 
4.62 


3.18 
4.01 


2.76 
2.61 


2.29 
1.96 


1.93 
1.68 


1.60 
1.47 


1.21 
1.10 


1.13 
1.27 


1.39 
1.38 


1.64 
1.60 


1.96 
1.68 


2.48 
2.28 


3.12 
2.52 


3.66 
2.84 


4.24 
3.39 


Exp. 
Obs. 


ee 
ee 


16.07 
19.37 


11.53 
13.93 


8.08 
10.08 


6.30 
9.12 


4.84 
6.31 


4.05 
6.65 


3.60 
4.78 


3.12 
4.95 


2.65 
4.52 


1.92 
2.61 


1.82 
2.88 


2.11 
3.35 


2.63 
3.16 


3.07 
4.82 


3.66 
4.55 


4.76 
7.12 


5.84 
8.26 


6.52 
9.24 


Exp. 
Obs. 


nil 


15.09 
25.84 


9.97 
18.83 


7.08 
9.93 


5.44 
8.34 


4.56 
7.01 


3.92 
7.11 


3.37 
5.37 


2.93 
4.45 


2.60 
3.88 


1.99 
2.99 


1.83 
2.31 


2.29 
3.22 


2.72 
3.79 


3.24 
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4.16 
5.16 


5.08 
5.78 


5.68 
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8.42 
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FIG. 16: [color online] The observed LLR (solid line) as a function of Mr using (a) all channels, (b) e/x channel, (c) ee 
channel, (d) nji channel. Also shown are the expected LLR distributions for the background only hypothesis (dashed line) and 
for the signal+background (dash-dotted line) hypothesis, with the dark grey and light grey bands indicating ±1 and ±2 s.d. 
fluctuations of the expected LLR distributions for the background-only hypothesis, respectively. 
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